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We show that when the standard techniques for calculating fractal dimensions 
in empirical data (such as the box counting) are applied on uniformly random 
structures, apparent fractal behavior is observed in a range between physically 
relevant cutoffs. This range, spanning between one and two decades for densities of 
0.1 and lower, is in good agreement with the typical range observed in experiments. 
The dimensions are not universal and depend on density. Our observations are 
applicable to spatial, temporal and spectral random structures, all with non-zero 
measure. Fat fractal analysis does not seem to add information over routine fractal 
analysis procedures. Most significantly, we find that this apparent fractal behavior 
is robust even to the presence of moderate correlations. We thus propose that 
apparent fractal behavior observed experimentally over a limited range in some 
systems, may often have its origin in underlying randomness. 



1 Introduction 

Fractal structures have been observed ip,-a large variety of experimental sys- 
tems in physics, chemistry and biologyLH3 Unlike exact (mathematical) frac- 
tals which are constructed to maintain scale invariance over many orders of 
magnitude, and most existing physical models displaying fractal behavior0 for 
empirical fractals the range over which they obey a scaling law is necessarily 
restricted by upper and lower cutoffs. In most experimental situations this 
range may be quite small, namely not more than one or two orders of magni- 
tude (FigJlh. Nevertheless, even in these cases the fractal analysis condenses 
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Figure 1: To obtain a general idea about the experimental status of fractal dimension 
measurements we collected all such measurements presented in Ref.[2] and measured the 
width of the linear range in the log-log plots (measured in decades) over which the FD was 
determined. This histogram shows the number of plots as a function of the number of decades 
of the linear range. One can see that most experimental measurements of fractal dimensions 
are based on data that extends between one and two decades. Note further that all the data 
with three and four decades, which come from a single paper, is the determination of the 
Hurst exponent for temporal and not structural data. 

data into useful relations between different quantities and often provides valu- 
able insight. 

Motivated by the yet inexplicable abundance of reported fractals, we con- 
sider here the apparent fractal properties of systems which are governed by 
uniformly random distributions. The reasons for this choice are several. First, 
randomness is abundant in nature. Second, although a uniformly random sys- 
tem cannot be fully scale invariant, it may, as we show below, display apparent 
fractality over a limited range, perhaps in better agreement with the actual 
ranges observed than a model which is inherently scale free. Third, a model of 
uniform randomness is a convenient limit, on top of which correlations can be 
introduced as perturbations. 

2 The Basic Model 

To illustrate our ideas we use a model that consists of a random distribution 
of spheres of diameter d, in the limit of low volume fraction occupied by the 
spheres. The positions of the centers of these spheres are determined by a 
uniform random distribution and the spheres are allowed to overlap. This 
model may approximately describe the spatial distribution of objects such as 
pores in porous media, craters on the moon, droplets in a cloud and adsorbates 
on a substrate as well as some energy spectra and random temporal signals. 

To simplify the analysis we consider here (without loss of generality) the 
one dimensional case, where the spheres are M rods of length d which are 



2 



11=0.01 , log, (d/L)=-6 




-8.0 -7.0 -6.0 -5.0 -4.0 -3.0 -2.0 
l°g,„(r) 

Figure 2: Comparison of simulation results (circles) to the theoretical prediction of Eq.(^) 
(solid line) for the number of intersected boxes as a function of their size, for one dimensional 
penetrable rods. The coverage is r\ = 0.01 and the rod length is d/L = 10" 6 . The cutoffs 
are manifested as the two knees in the graph. The lower bound ro is seen to be located at 
r = d/L. The upper bound r\ is at r = (l/»? — l)d/L, also conforming with the prediction 
in the text. Also indicated is the estimated middle point r e . 



placed on a line of length L 3> d. The positions of the rod centers are deter- 
mined by a uniform random distribution. The rods are allowed to overlap and 
are positioned with no correlations. An information-theory argument can be 
used to show that this distribution is generic, or "minimal" , in the sense that 
it is characteristic of physical processes in which only the first moment (such 
as the density) is determined from outsidetm Below we calculate the fractal 
dimension (FD) of the resulting set using the box-counting (BC) procedure, 
which is a common techniques for determination of FD in empirical dataE2l In 
the BC technique one divides the embedding space into boxes of linear size I. 
It is convenient to work with the dimensionless quantity r = l/L for the box 
size. The number of boxes that have intersection with the measured object, 
N(r), is then plotted vs. r on a log- log scale. The range of r is limited from 
below by the finest feature in the object and from above by the entire object 
size. Apparent fractal behavior is commonly declared in a range bound be- 
tween physical cutoffs if the log-log plot of N(r) vs. r is linear over one or 
more decadesLil in that range. The dimension is given by: 

D = -slope {log(r),log[iV(r)]}. (1) 

We will now show that our model generates approximate linearity over a range 
which would conventionally be accepted to indicate fractality. The lower cutoff 
is given by the rod length, 

ro = d/L, (2) 
since below this scale no new information is obtained by decreasing the box 
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Figure 3: Simulation results (circles) for the number of intersected boxes N(r) vs. r in 
the experimentally relevant range (between cutoffs), along with a linear regression fit for 
coverage r) = 0.1 (d/L = fCT 5 , M = 10 4 ). 

size. The upper cutoff is determined by the average distance between adjacent 
rod edges, 

n = l/M-d/L, (3) 

because above this scale (on averages-all boxes are occupied. This allows us to 
define an estimated scaling range asa 

A e = log 10 (ri) -log 10 (>o). (4) 

Since the actual value of N(r) depends on the particular set of random numbers 
drawn, one can only obtain the expectation value (N(r)). However, the law of 
large numbers ensures that for a large enough number of rods, the deviations 
from this value will be insignificant. _ 

Following probabilistic arguments of the type used by Weissbergp and 
Torquato and Stell[L3 one obtains_Lhat out of the total of 1 jr boxes the number 
of boxes that intersect the set isEJ 

(N(r)) = ±{l-{l-(r + d/L)} M }. (5) 

Simulation results in terms of the coverage rj = Md/L are shown in Fig.^j, along 
with the theoretical prediction of Eq.(||). An excellent agreement is evidcnttJ 
Next, we examine the apparent FD [Eq. ([[])], by mimicing the standard 
experimental procedure of using linear regression analysis between the cutoffs. 
The simulation results and the linear fit for r\ = 0.1 are shown in Fig.|^ for 
the range which is used to determine empirical FDs. More than a decade 
of linearity is observed for this high coverage. The slight inflexion of the 
simulation results may be smeared out by noise in a real experiment. We next 
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Figure 4: The range of linearity, A, as a function of imposed coefficient of determination, 
R 2 , in a linear regression analysis. 



evaluate the slopes and actual ranges of linearity A (generally ^ A e ), under 
varying degrees of strictness of linearity, as measured by the coefficient of 
determination R 2 . Typical results are shown in Fig.^J, where, e.g. for 7] = 0.01, 
more than two decades of linear behavior are exhibited for a required value 
of R 2 of below 0.975. This is well within the experimental norm as most 
experimental measurements of fractal objects do not extend for more than two 
orders of magnitude (Fig.|l|). Moreover, this agreement with experimental data 
is in contrast to that of most other physical models of fractality, which predict 
much larger rangesQ Increasing r/ beyond 0.1 results in a decline of both A 
and A e to below one decade and hence the apparent fractality is restricted to 
i] < 0.1. 

The results of the regression analysis for the apparent FD as a function 
of rj are shown in Fig.|| and are further compared to an analytical expression, 
obtained by calculating the logarithmic derivative of N(r) at the estimated 
middle point r e = a/totT, in the M — > 00, constant coverage limitta 

D = l . (6) 

exp (77 + 07(1 - 7])j - 1 

As seen in Fig jHJ, the FD predicted by Eq.(|^) is somewhat lower than the 
regression result and can serve as a lower bound. In the limit of small 77, one 
can further simplify Eq.(|^) and obtain 
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Figure 5: Apparent fractality (FD) as computed by linearregressipn with R 2 = 0.995 (upper 
curve). The predictions of the analytical equations, Eqs.ffej) and (0) (for the penetrable and 
impenetrable rods) are accurate lower bounds and differ only marginally (two overlapping 
lower curves). This indicates the dominance of randomness over correlations. 



3 Impenetrable Rods 

To examine the effect of correlations on the apparent FD we next consider a 
model in which rods are randomly located as before but with the restriction 
that the rods cannot overlap. The system is assumed to be at equilibrium. 
The excluded volume effect clearly creates correlations in the positions of the 
rods. This example is also fully solvable! and represents an important class 
of systems with correlations such as models of hard-sphere liquids and energy 
spectra with level repulsion. We will now show that the correlation introduced 
by the non-overlap restriction merely modifies the apparent fractal character 
of the system. For this case, the expected number of intersected boxes isH 

(N(r)) = ± (l - (1 - „) (l - . (8) 

Fig.|| shows the number of intersected boxes (N(r)) vs. r both with [Eq.mh] 
and without [Eq.(||)] overlap. The behavior in the two cases is qualitatively 
similar and virtually indistinguishable for low coverages. Fig thus demon- 
strates that the apparent fractal behavior due to randomness is only slightly 
modified by moderate correlations. As in the overlapping rods case, we can 
now use Eq.(Q) (with the slope calculated at r — r e ) to calculate a lower bound 
for the apparent FD. The result (for large M), 



D = l . V^l^ 1 (9 ) 



exp 



is shown in Fig.^. The important observation is that for a broad range of 
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Figure 6: Comparison of box-counting predictions in penetrable and impenetrable rods 
cases. The results for penetrable [Eq. (|j)] and impenetrable rods [Eq.(jfj)] virtually coincide 
for rj < 10 — 2 (lower two curves). For rj = 0.1 a barely noticeable difference develops (upper 
two curves). In both cases d/L = 10 — 6 . 



low coverages the apparent FD's of penetrable and impenetrable rods nearly 
overlap. This is the relevant range for fractal measurements and therefore 
we find that correlations of the type considered here have little effect on the 
apparent fractal nature of the system. 

4 Fat-Fractal Analysis 

In this section we treat the penetrable spheres model for the case of two- 
dimensional (2D) disks, from the point of view of fat-fractal analysis. A 
fat fractalis defined as "A set with a fractal boundary and finite Lebesgue 
measure" £3 The fat-fractal approach is natural for our model, since the set 
of disks clearly has non-zero measure. Fat-fractal analysis can be performed 
on experimental data (but rarely is) in those cases where the resolution of 
the measurement device is finer than the lower cut-off, which is required for.-a. 
knowledge of the measure of the studied set. An example is helium scatteringlia 
In the present case we show that the measure of the set of disks can be found 
analytically. In order to measure the fat-fractal scaling exponent 7, one per- 
forms, as in the standard fractal analysis, a box-counting procedure: 

= Hm log[A(r)} ^ r2N{r) _ ^ 

r^r' lOg(r) 

where fi is the normalized Lebesgue measure of the set. The fractal dimension 
itself is given by 

D ff = 2- r (11) 
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Figure 7: Fat fractal analysis of random disks. Analytical results of Eq. Ji"o|) for A(r) are 
shown at three combinations of coverages^ and disk numbers M . Circles indicate the 
positions of the cutoffs according to Eqs.(H),(H). Inset: Linear regression coefficient 7Z 2 for 
regression in-between the cutoffs. 

In the nonlinear dynamical systems literature, where fat fractals were first 
introduced pi r* — 0. In the context of real-space sets, there exists a lower 
cutoff ro > 0, and hence also r* > 0. One should bear this in mind whenever 
fractal theory is applied to real-space systems with an inherent non-vanishing 
smallest scale. 

Consider then again a system of M uniformly randomly positioned disks 
of equal radius R — d/2, located at low 2D coverage 772 given by: 

r/2 = MnR 2 /L 2 = (7r/4)77?, (12) 
on a surface of area L 2 . The effective "ID coverage" 

77! = VM2R/L, (13) 

is defined for convenience of comparison with results in ID and 3D. In order 
to find /Lto, imagine that the surface is initially empty, and randomly choose a 
point on it. Next locate a disk of radius R at a random position on the surface. 
The probability that it does not include the chosen point is proportional to 
the free area, namely q\ = (L 2 — irR 2 )/L 2 . The next disk is also positioned 
completely randomly, so that the probability for the point to be outside of 
both disks is just q\. Clearly, after random placement of M disks, the point 
will lie in the uncovered region with probability qf 4 , and therefore will be in 
the disk-covered region with probability pm = 1 — (1 — ttR 2 / L 2 ) M . On the 
other hand, this probability is just the expectat ion v alue of the normalized 
disk union area, /j,q/L 2 . Thus for large enough MEjEj 



Mo = [1 - (1 - nR 2 /L 2 ) M ] L 2 . 



(14) 
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Figure 8: Analytical and regression slope between the cutoffs (dashed and long dashed) of 
fat fractal analysis, and 2 — D of "thin" fractal analysis. 



A modified argument can be used to evaluate the BC function for our basic 
modeltl The result for the expected number of occupied boxes is: 



N(r) 



L' 2 



1 - (1 - r 2 - irR/L - Tr{R/Lf 



, AT 



(15) 



Simulation:!! (not shown here) confirm the ID version of this result to excellent 
accuracy for M as small as 100. Taken together, Eqs.(|T^),(^5|) determine the 
fat-fractal exponent 7, using Eq.(10). Analytical results are shown in Fig.^for 
three rji/M pairs. The effect of changing M at constant coverage (solid and 
short-dashed lines) is a rigid translation of the curve in the plane. This implies 
that the coverage is the important parameter in determining the slope i.e., 
the FD. Circles indicate the positions of the cutoffs according to Eqs.(g),@. 
Beyond the lower cutoff the slope tends to 1, beyond the upper cutoff - to 0. 
In-between the cutoffs, a nearly straight line is observed, in agreement with 
apparent fractal behavior. In order to find 7 it remains to determine the point 
r*. For disks, in analogy to the discussion for rods, the cutoffs are given by: 



r = 2R/L 

n = 1/VM -2R/L. 



(16) 



As in Sec.|| we choose r* as the estimated middle point of the scaling range, 



r e = \/rori, 



(17) 



and find 7 by evaluating the logarithmic derivative of A(r) there. The result 
is: 
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Figure 9: Log-log plot of A(r) = r n N(r) — fio (solid line) and just r n N(r) (dashed line). 
Inside the physical range, the measure fig has little effect. Circles indicate the cutoffs. 



d\og[A(r)} 27?i (i - m + Vm-vi) (18) 

dlog(r) re exp[77i(l- jji+2071- 77?)] -l' 

This result is compared in F-jg.|| to the result of the regular ("thin") fractal 
analysis, for which we foundB 



2 D= ^[7710771 - 77? + 2771 -rj 2 } 

cxp{77 2 + \ [771 077i - 77? + I(277i - 77?)]} - 1 ' 

The two curves differ only slightly. The analytical fat-fractal result is also 
compared in Fig.|| to the procedure followed in typical experimental analysis 
of fractal scaling data: a linear regression between the physical cutoffs (tq and 
r\ in our case). The trend is similar, and the agreement is quite good for the 
higher coverages. In any case the analytical Eq.(|l8|) serves as an accurate upper 
bound to the expected regression result for 7. The corresponding regression 
coefficient 1Z 2 (inset of Fig.fj]) does not fall below 0.9985 which indicates a very 
high quality regression, certainly by experimental standards. Note that 1Z 2 
remains very high even for 771 < 10~ 2 (i.e., a scaling range A e > 2). This is 
a wider range than found for the ID version of "thin" fractal analysis. There 
the apparent fractality was observed in a range of 1-2 decades if 77 < 10 _1 
and 1Z 2 > 0.97 are required (Fig.^). This range improvement is not a direct 
outcome of the inherent fractality we found between the cutoffs, but is due to 



the differences in the response of the linear regression procedure to Eq.(10), in 
comparison to Eq. ([[]): It is the multiplication of the BC function by r 2 in the 
former which is responsible for the improved range effect, through the increase 
in the slope of the log- log plot. One should be cautious therefore to eliminate 
slope-biases in analyses of scaling properties on a log-log plots. Furthermore, 
as clearly seen in Fig.|l the essence of the fat fractal analysis, namely the 
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subtraction of the measure no of the set, has essentially no effect on the slope 
and on the location of the cutoffs and thus does not provide us in this case with 
added information. Being left then with the choice between N(r) and r 2 iV(r), 
there does not seem to be a clear reason to opt for the latter. We conclude 
that for low density systems, such as in this report, fat fractal analysis is not 
necessary. 

5 Conclusions 

In summary, we have shown that random structures, which are generic in 
experimental situations where only the first moment of a distribution is deter- 
mined, give rise to apparent fractal behavior within physically relevant cutoffs, 
with a non- universal FD. Although this is not a mathematically rigorous frac- 
tality, in the sense that the scaling is not strictly a power law, it is a physical 
fractality: It satisfies the conditions of high-quality linear regression in the 
physically relevant range of observation. Since experiments rarely observe a 
perfect power law, we believe that the possibility of approximate scaling should 
be considered in theoretical models, if a more complete understanding of the 
experimental fractal data is to be achieved. It is likely that some of this data 
does in fact not reflect the existence of an exact power law, but rather an ap- 
proximate power law between cutoffs with a weak inflexion point in the log-log 
plot. The present model and its approximate scaling properties hint that this 
may be the case, e.g., for porous media. Moderate correlations have little effect 
on the apparent fractal properties and even in their presence it is still the un- 
derlying randomness that is the main contributor to the apparent power-law 
scaling relation. Elsewhere we showed that these results remain practically 
unchanged for higher dimensions and for a variety of size distribution profiles 
of the elementary building blocks!! We thus propose to consider randomness 
as a possible common source for apparent fractality. 
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